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Aim 

To understand the effect of indistinguishable particles on the 
statistics, and to apply this to the ideal gas. 

Objectives 

1. To explain the counting of states and derive the density of 
states. 

2. To explain the conditions under which the ideal gas model 
would not be valid. 

3. To derive the Maxwell-Boltzmann distribution. 

4. To derive the heat capacity for mono and diatomic gases. 

5. To explain the effects of ortho and parahydrogen. 
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4.1 Indistinguishable particles 
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Indistinguishable particles 



Indistinguishable particles are identical particles. 

It is not possible to tell one from the other. 

There is a special case. If they fixed in a solid, as with the 
paramagnet that we have seen previously, then we can tell 
them apart with the help of their location. In that case, 
however, they become distinguishable particles. 

We are therefore interested in cases when they are not fixed. 
Examples are atoms in a gas, electrons in a metal, and photons 
in a cavity. 

We shall start with the monatomic ideal gas. Each particle in 
this gas is a single atom. 
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How are they different from distinguishable particles? 



Suppose that there are 3 particles, each with 3 energy levels ei, 
62 and 63. 

Consider the macrostate ni = 1, = 1 and 77-3 = 1. I.e. each 
level has one atom: 




If the particles are distinguishable, then there are 6 microstates. 
However, if the particles are indistinguishable, there is only 1. 
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Just erase the numbers - they all look the same! 
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Ideal gas 



There are a few steps to deriving the statistics for an ideal gas: 

1. We start with a single atom in a 3-D box. We solve 
Schrodinger's equation and obtain the energy states. 

2. At room temperature, these energy states would be very 
close together. They get even closer at higher energy. We need 
to find the density of states, which gives the number states per 
unit energy. 

3. It turns out that the large density of energy states means 
that it is very unlikely for any two atoms to even occupy the 
same levels. This makes it possible to simplify the calculation 
and obtain the energy distribution of the atoms. 
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4.2 Counting states 
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A particle in a 3-D box 



We first determine the energy states for a particle in a 3-D box. 





We solve Schrodinger's equation for the particle: 

2m \dx 2 dy 2 dz 2 J 

where V is zero in the box, and infinite at the walls. So the 
wave function must have zero amplitude there. 
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Energy states of the particle 



This is done in the standard way by trying the solution 

if; = sin k x xs\n k y y sin k z z 

where (k x ,k y ,k z ) are unknown constants. They can be 
determined from the boundary condition that the wave function 
must have zero amplitude at the wall. 

In the x direction, suppose the walls are at x = 0 and x = a. 
The condition can be written as: 

sin k x a = 0. 

This gives 

7 7 

k x a = n x 7r or k x = , 

a 

where n x is a positive integer. This means that k x is quantised. 
Similarly, for the y and the z directions, we find 

k y = and k z = . 

a a 
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Energy states of the particle 



We can then substitute the solution 

ip = si n k x x si n k y y si n k z z 
into the Schrodinger's equation 

2m \dx 2 dy 2 dz 2 J 
to find the energy. This gives 

H 2 2 2 2 
— {k x + k y + k z ) = E. 

It can also be written in the form 

h 2 k 2 

= E 



2m 

where 



k — ^ z ' 



k is called the wavevector. 
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Energy states of the particle 



We have obtained the energy of the particle: 

,2 , 7 2 , 7 2 



h2 

E — — {k x + k y -\- k z ) 



2m 

and the components of the wavevector: 

_ n x iT _ n y 7r _ n z 7r 

Kx 5 rvy 3 nO Kz — • 

a a a 

Substituting, we get 

h 2 

2 2 2 
£n x n y n z = — o( n x n y ~l~ n z)' 

E is replaced by the symbol e nx n y n z - This is to indicate that it 
is quantised, and depends on the positive integers (n x ,n y ,n z ). 
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Counting states in 2-D 



We need a way to determine the number of states below a 
certain energy, s. This would then allow us to find the density 
of states. 



To develop the idea, start from the simpler case in 2-D. If we 
consider a particle in a 2-D box, the energy would be 

h 2 2 2 

£ n x n y — - o( n x + n y) 

Notice that n z has been left out. We can represent the possible 
energy states by treating (n x ,n y ) as co-ordinates, and plotting 
them on a graph: 
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Counting states in 2-D 




Remember that (n x ,n y ) must be positive integers. So the 
points are 1 unit distance apart. 

If we draw a circle of radius R, it would be described by the 
equation 
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Counting states in 2-D 
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Compare the equations for the circle and for the energy: 

ix >2 _ _2^_2 

8ma 2 e n , 



R = n x + n y 



^2 — n x\ n y 



If we associate the energy with the radius: R 2 = 8ma z s/h 



we could find the number of states below that s, by counting 
the number of points enclosed by the circle. 
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Counting states in 2-D 



• 4 



—I 1 1 1— 



Note that all the points are in the first quadrant only. 

Note also we shall assume R is much larger than the distance 
between adjacent points. (We would need to justify this later.) 

The distance between adjacent points is 1 unit. Consider a unit 
square centred at the point. 
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Counting states in 2-D 



If we draw such a square for every point, the squares would fill 
up the whole quadrant. 
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Instead of counting the points, we can count the squares. 



Since the area of each square is 1, the number of squares 
would just be the area of the first quadrant, 7rR 2 /4. 
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Counting states in 2-D 



—I 1 1 1 1 1 L| 1 1 1_ 



IR "* 



Since the number of states is the number of points, this 
number is also 7tR 2 /4. 



Recall the relation between energy and radius: R 2 = 8ma 2 s/h 2 



Therefore, in 2-D, the number of states below e is: 

= 27rma s/h . 
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Counting states in 3-D 



In 3-D, the energy states are given by: 

Sma s n 



x n y _ 2 I 2 i 2 
= n% + riy + 



Compare this with the equation of a sphere. 




Since (n x ,n y ,n z ) are positive, the states are given by the points 
in the first octant (1/8) of the sphere. 
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Counting states in 3-D 



Next, we need to work out how to count the states in 3-D. We 
need to make the following changes: 

1. (n x ,n y ) to (n x ,n y ,nz) 

2. Points in 2-D to points in 3-D 

3. Circle to sphere. 

4. Quadrant to octant (1/8 of a sphere). 

5. Unit square to unit cube. 

6. Area to volume. 
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Counting states in 3-D 



The formulae are as follows: 



1. 1/8 of the volume of a sphere is | x ^nR 3 . 

2. The relation between energy and radius remains the same: 

2 _ 8ma 2 e 

R 

3. In 3-D, a unit cube is associated with each state. So the 
number of states is the volume of the octant: 

, ^ 14 f8ma 2 s x3/2 
G(e) = - x -7r o — 

w 8 3 V ^ 
= 3^ (2m£) 7 

where V = a 3 . 
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Density of states 



We can now find the density of states. We have found that the 
total number below energy e is 

If we increased the energy by 5s, this number would change by 
SG(s). 

The density of state is the number of states per unit energy. 
Therefore this is given by 5G(s)/5s, or 

g (e) = dG(s)/ds = 4 ^(2m£) 1 / 2 
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4.3 Ideal gas 
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Ideal gas energy states 



We should now justify our assumption that it is very unlikely for 
two atoms to occupy the same energy state. 



At a temperature T, we know from simple kinetic theory of 
ideal gas that the average energy of a gas atom is about 
3k B T/2. 

First, we need to know how many energy states there are 
below this energy at room temperature. We can use the 
formula we have derived: 



G(e) = f^(2me)3/2 



by setting s = 3k B T/2 
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Ideal gas energy states 



G(3k B T/2) = ^(2me) 3 / 2 

Putting in the numbers for one mole of ideal gas at room 
temperature, we find that the number of states below 3k B T/2 
is about 10 30 . 

1 mole of gas contains about 10 24 atoms. 

This means we have about 10 30 10 24 = 10 6 states per atom 
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Ideal gas energy states 



This means we have about 1 million energy states for every 
atom. 



CD 
i_ 

<D 

C 

LU 



So it is extremely unlikely that two atoms would ever occupy 
the same energy state. 

We are now going to derive the energy distribution of these 
atoms. 
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Grouping energy levels into bundles 



We shall borrow the idea from histograms - divide the energy 
levels into bins, or bundles. 



D) 
O 

c 

LU 



g j+1 levels in the (/+1)" 7 bin 
g, levels in the I th bin 
g t i levels in the (/-1) th bin 



Suppose each bundle covers an energy interval ds. 

We want to find out how many ways we can arrange a number 
of atoms in a particular bundle. 
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Energy levels as partitions 



Consider the i th bundle. There are gi energy states and 
atoms. 



C 
LU 



Recall that it is extremely unlikely for two atoms to be in the 
same state. 



Think of each state as a partition. The partition and the atom 
are like two different objects. We need to find the number of 
ways to arrange ni atoms and gi partitions in a row. 
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Number of arrangements 



There are atoms and gi partitions. 



We can use the combination formula n C r = 



This gives, for one bundle, ( n *~j~^) ! ways 



r\(n—r) ! 



To find the total number of arrangements for all bundles, we 
must multiply together the answer for every bundle: 

Q = n (n '|^ )! C 1 ) 



Since the number of energy states in each bundle is much 
larger than the number of atoms, gi » r^, we can simplify this. 
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Expand 



q _ jj (gz + nj)...(gi + l)(flfj...(l) _ jj + ni)...{gi + 1) 



Since ^ >> n if + « ^ and + 1) « ^. Therefore 



This gives the number of microstates in the macrostate n. 
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If a macrostate is specified by (711,712,713,...), or ni in short, 
then the above expression gives the number of microstates in 
the macrostate. 

We are interested in the most likely distribution of 
(711,712,713,...), or the most probable macrostate. 

We can once again apply the Lagrange multiplier to find this. 
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Most probable macrostate for the ideal gas 



We apply the Lagrange multiplier method as follows: 

1. Instead of Q, maximise InQ. 

2. Apply Stirling's approximation: 

In Q = Y,Wi ln 9i ~ ™i ln ™i + rii] = n i 'n(^M)] + N 

i i 

3. The number of particles is fixed: 

i 

4. The total energy is fixed: 

i 

5. Define the Lagrange function as In Q + XiN + X 2 U. 

6. To find the maximum, differentiate with respect to nf 

In gi - In ni + Ai + X 2 £i = 0 
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Most probable macrostate for the ideal gas 



We have found 

In g i - In rii + Ai + X 2 £i = 0 
Rearranging, we get 

rii = A^exp(A 2 ^), 

where 

A = exp(Ai). 

In the same way as before, by using our knowledge of entropy, 
we can show that 

1 



A 2 = - 

k B T 

This gives the Boltzmann distribution again: 

rii = Agi exp(-Si/k B T), 
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Approximate to continuous energy 



We are mainly interested in the energy range of the order of 
k B T (which means up to a few times of that). Most of the 
atoms would fall in this range. 

We have seen previously that there is a very large number of 
energy states in this range. 

Since the energy levels are so closely spaced, we can 
approximate them to continuous energy. 

We have previously obtained the density of states, which is the 
number of states per unit energy, g(s). 

Consider the i th bundle. Assume that this occupies an energy 
interval of de. 

The number of energy states is then given by gi = g(s)ds. 
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Approximate to continuous energy 

In a similar way, we can define a number density n(s) as the 
number of particles per unit energy. 

So the number of atoms in the bundle would be 

Therefore, the distribution that we have obtained, 

rii = Ag i exp(-f3e i ), 

can be rewritten as 

n(s)ds = Ag(e) exp(—Si/k B T)ds 
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Approximate to continuous energy 



Substituting the formula for g(s), we get 

, N A A-mirV , ,1/9 
n(s) = A—^&ms) 1 / 2 exp(-Si/k B T) 

Note that the total number of atoms is then given by 
integrating over all energy: 

•oo 



roo 

N = / n(s)ds 
JO 



This relation is known as the normalisation condition 



PHYS393/Hock 



36 



2009-2010 



Maxwell-Boltzmann distribution 



The unknown constant A can be determined by integrating the 
normalisation condition (with the help of the table of integrals): 

A= N( h 2 \ 3/2 
~ V \2irmk B TJ 



This gives 



»W = (^372 W 1/2 -p(-, A B T) 



The energy is just the kinetic energy, e = mv 2 /2. So 

2\V2 / m \ 3/2 
k B T t 

This is the Maxwell-Boltzmann distribution of speeds. 



/2\V2 / m \ 3/2 o 
n(v) = N (-J I — J exp(-rm; 2 /2£^TV 
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In case you are interested 



There is actually a simpler way to derive the 
Maxwell-Boltzmann distribution, without using Lagrange 
multipliers. 

You can find it here: 

http : / / galileo . phys . Virginia . edu/ classes/ 252/kinet ic_theory . html 
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At the Columbia University in New York, 



Miller and Kusch measured the speeds of the atoms in 
potassium vapour. They agree very well with the formula 
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Fig. 2. Schematic diagram of the apparatus designed 
to measure velocity distributions. 
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Fig. 4. Typical potassium velocity distributions. The vapor 
pressures in tbe ovens are given in Table IT. Run 31 was made 
with thick oven slits and runs 57 and 60 with thin slits. 



This was published in Physical Review in 1955. 
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When would it break down? 



We have assumed that: "It is extremely unlikely that two 
atoms would ever occupy the same energy state." 



O) 

C 

LU 



Recall that we have arrived at this by looking at the number of 
states below the mean energy 3k B T/2, when T is at room 
temperature. 



If T is very small, there may be far fewer states below the mean 
energy. Then it would be very likely for two atoms to occupy 
the same state, and the Maxwell-Boltzmann distribution would 
not be valid. 



Instead, quantum statistics - such as Fermi-Dirac or 
Bose-Einstein - have to be used. We shall learn about these 
later. 



PHYS393/Hock 



40 



2009-2010 



A note on the recommended text 



If you are using the recommended text: 

Statistical Mechanics - A Survival Guide, 
by A. M. Glazer and J. S. Wark 

please note that the Lagrange multiplier method given on page 
49 is incorrect. 

Please use the version given in these lectures. 
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4.4 Partition function 
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The partition function 



The partition function would be useful for deriving the 
macroscopic quantities like heat capacity, pressure and entropy 
from the microscopic properties. 

It would also be useful for extending the calculation to diatomic 
gases, and include effects like rotation, vibration and even 
nuclear spin states. 

In the lectures on paramagnetic salts, we have derived formulae 
for energy and entropy in terms of the partition function. The 
same formulae, unfortunately, cannot be used for the ideal gas. 
The reason is related to the fact that the particles in the ideal 
gas are indistinguishable. 
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The partition function 



To see how to to obtain the partition function for the ideal gas, 
we first rewrite the formulae for the paramagnetic salts: 

Total energy 

9 <9lnZ 

U = NkuT 

B 8T 

Entropy 

d\nZ 



S = Nk B \nZ + Nk B T- 

ol 

Helmholtz free energy 

F = -Nk B T\nZ 

Notice all of these formulae contain the term AMnZ. Rewrite 
this as 



N\nZ = \nZ N 



Z N is called the system partition function. Denote this by Z N 
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The partition function 



In terms of the system partition function, the formulae can be 
written as follows: 



Total energy 



o<9ln Z N 
U = k B T 2 



dT 
Entropy 

S = k B \nZ N + k B T 9U ^ 
Helmholtz free energy 

F = -k B T\nZ N 

In the case of the ideal gas, the system partition function has 
to be replaced by 

_ Z N 
Zat — 



' N N\ 

The reason is related to the fact that the particles are 
indistinguishable. Dividing by 7V! would avoid overcounting. 
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The partition function 



In terms of the system partition function, we can continue use 
the same formulae for energy, entropy, and free energy. 

The difference is that, for distinguishable particles, the system 
partition function is 

Z N — Z N , 

whereas for indistinguishable particles, the system partition 
function is 

_ Z N 

Zat 



' N N\ 

For convenience, we may refer to the original partition function 
Z as the single particle partition function. 



PHYS393/Hock 



46 



2009-2010 



The partition function 



We have previously written the single particle partition function 
as a summation: 



Z = $>xp - 



For the ideal gas, the spacing between energy levels are small 
compared to k B T. So the summation may be approximated by 
an integration: 



Z — I o(s)exp( e —\de 

^ V k B Tj 
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The heat capacity 



Substituting the density of states, 

00 AmnV 



Z = 



(2ms) 1 / 2 exp ( — ) de 

V k B Tj 



r 0 h 3 

This can be integrated with the help of the table of integrals: 

z _ v ( W ) 3 /2 

Recall that there is a set of formulae that we can use to find 
the energy and entropy from the partition function. The energy 
is then given by this formula: 

o<9ln Z N 
U = k B T 2 

Substituting the system partition function Z N /N\, we find 

v = 3Nk B T 
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The partition function for indistinguishable particles 



The heat capacity is then obtained by differentiating with 
respect to temperature: 

c = ™ = 3 -N kB . 

dT 2 a 

This agrees with the result from kinetic theory of the ideal gas, 

http : //en . wikipedia . org/wiki/Kinetic_theory . 
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Ideal gas law 



Next, we shall show that we can derive the ideal gas law 

P V = nRT 

using the partition function. In order to do so, we shall use the 
Helmholtz free energy to calculate the pressure p. 

Recall from the lectures on paramagnetic salts that the 
differential of the Helmholtz free energy is: 

dF = -SdT - pdV. 
which leads to this formula for p\ 
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Ideal gas law 



Here is a list of all the formulae we would need: 

To find the pressure from the Helmholtz free energy: 

(-) =-p- 

\dV/T 

The partition function for the Helmholtz free energy: 

F = -k B T\nZ N , 
The system partition function for indistinguishable particles: 

_ Z N 



' N N\ 

The partition function for the ideal gas: 

z = v 
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Ideal gas law 



Substituting the system partition function 

_ Z N 

into the Helmholtz free energy formula 

F = -k B T\nZ N , 

we find 

F = -k B T(N\nZ - In N\). 
To see what InZ looks like, we substitute the partition function 




This gives: 



In Z = In V + terms not containing V. 
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Ideal gas law 



In order to find the pressure p, we need to take the partial 
derivative with respect to V: 



(-) = 

wvJt 



■p- 



So terms that do not contain V would vanish. Substituting the 
expression for InZ into the equation for the Helmholtz free 
energy formula, we find 



F = —Nk B T\n V + terms not containing V. 

Differentiating with respect to V, we get 

( dF \ _ _Nk B T 
\dVJT ~ 

Since this partial derivative is equal to — p, we have 



_ Nk B T 

Rearranging, we have the ideal gas law: 

P V = Nk B T 
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Ideal gas law 



To write this in the more conventional form, we know that the 
Boltzmann constant is 

k B = R 



Na 

where R is the gas constant and is the Avogadro constant. 
We also know that the number of moles TV is related to the 
number of particles n by 

TV = nN A . 

Substituting these, into the above equation for p, we get the 
more conventional form of the ideal gas law: 

P V = nRT. 
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4.5 Boltzmann constant 
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Boltzmann constant 



We can now return to the earlier lectures on derivation of the 
Boltzmann distribution 

rii = Aexp(si/k B T) 
and tie up a loose end. 

Recall that the kgT in the Boltzmann factor comes from one of 
the Lagrange multipliers: 

1 



A 2 = 

k B T 

We have seen how this is related to temperature T using our 
knowledge on entropy. 
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Boltzmann constant 



We have also seen that there should be a constant, 



k B T 

but we have not looked at why the constant should be the 
Boltzmann constant 

k B = R 



N A 

where R is the gas constant and N A is the Avogadro's 
constant. We are ready to do that now. 

Recall that we have derived the following form of the ideal gas 
law: 

P V = Nk B T. 

We have derived this from the formula for the Helmholtz free 
energy expressed in terms of the partition function. So the 
Boltzmann constant k B here originates from the Lagrange 
multiplier. 

PHYS393/Hock 57 2009-2010 



Boltzmann constant 



The ideal gas law, on the other hand, may be viewed as an 
empirical law: 

P V = nRT. 

Under suitable conditions, the gas constant R can be measured 
using real gases, and has the value of 8.314 J mol -1 K _1 . 

Think of kg as an unknown constant from the Lagrange 
multiplier. We can find it by comparing the empirical law above, 
with the gas law that is derived using the partition function: 

P V = Nk B T. 
Comparing the above 2 equations, we find 

Nk B = nR. 
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Boltzmann constant 



Since the number of moles n is related to the number of 
particles TV by 

TV 

n = , 

we have 

Nk B = —R 

N A 

and finally, we obtain the value for k B : 

t R 

k B = 



which is indeed the Boltzmann constant as we have assumed 



We should note that the ability to derive the ideal gas law from 
the partition function is one of the successes of statistical 
mechanics. It is evidence that our rather abstract ideas about 
probabilities, microstates and Lagrange multipliers do give a 
realistic description of the actual behaviour of real gases. 
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4.6 Diatomic ideal gas 
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Diatomic ideal gas 



The molecule of a diatomic gas has two atoms bonded 
together. Additional types of movements, or degrees of 
freedom, are possible. 



For the monatomic gas, the atom has only translation motion. 
In the diatomic gas, the molecule also has rotation and 
vibration. 

The total energy is: s tota | = s trans + ^rot + ^vib 

More details about diatomic molecules can be found here: 




translate 



rotate 



vibrate 



http : //en . wikipedia . org/wiki/Diatomic_molecule 
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Paritition function for diatomic gas 



The probability of being in a particular total energy level can be 
found by multiplying the probabilities in the separate levels: 

Stotal exp(-£ to tal As? 1 ) = 

<?trans exp(-s trans //c B T) xg rot exp(-s rot /k B T) x # vib exp(-s vib //c B T) 

dtrans is the same density of states that we have seen previously 
for the monatomic gas. 

g rot and # vib are the density of states for rotation and vibration 
respectively. 

This means that the partition function can be written as 

7 _ 7 yZ yZ _ ( 2nmk B T \ 3 / 2 

^ — ^trans x ^rot x ^vib — Zj — v i ^ I x Z/ ro t x z, V j b 
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Rotational energy levels 



We know from quantum mechanics that the rotational energy 
levels are 

h 2 

sj = -J{J+l) 

where J is an integer and I the moment of inertia of the 
molecule. 

We also know that the components in the z direction are 
quantised, with quantum numbers: mj = — J, J — 1, +J. 

So there are (2J+ 1) states associated with the same energy 
level J. We say that there is a degeneracy of (2J+ 1). 

For example, for the angular momentum J = 1, there would be 
(2J+ 1) = 3 states with this angular momentum. These 
correspond to the components mj = -1, 0, 1. So we say that 
the angular momentum J = 1 has a degeneracy of 3. 
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Rotational energy levels 



The rotational partition function is then given by 

oo 

Z rot = £ (2J+ 1) exp[—h 2 J(J + l)/2Ik B T] 

J=0 

We can understand this by looking at the term J = 1 again. 
The Boltzmann factor for J = 1 would be 

exp[-fi 2 J(J+ l)/2Ik B T] = exv(-h 2 /Ik B T). 

Since there are (2J+ 1) = 3 states for the same J = 1, there 
are 3 states with the same Boltzmann factor. So the J = 1 
would contribute 3 times of the Boltzmann factor 
exp(—h 2 /Ik B T) to the partition function. 

In general, for the angular momentum J, there would be 
(2J+ 1) with the same Boltzmann factor. So this would 
contribute (2J+ 1) times of the Boltzmann factor 
exp[— h 2 J(J + l)/2Ik B T] to the partition function. 
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Excitation of rotation 



Since the rotation is quantised, it is quite possible that the 
molecule does not rotate at all if the is not enough energy. 

It would be useful to know if a molecule like hydrogen or 
oxygen would rotate at room temperature. 

The average rotation energy at a certain temperature is about 
k B T. How does this compare with the energy levels? 

The difference in energy between the J = 0 and J = 1 levels is 




I 

To have a reasonable fraction of the molecules in the first 
rotational state, this should be close to k B T. 
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The rotational temperature 



Equating As and k B T, we find that T is given by 

This called the "characteristic temperature of rotation," or the 
"rotational temperature." If the temperature is a lot lower than 
this, the molecules would hardly rotate. 

The moment of inertia for diatomic molecules may be 
estimated using / = mr 2 , where m is the mass of each atom 
and 2r the bond length. E.g. 



For nitrogen, 6 rot = 2.8 K 
For hydrogen, 0 rot = 82 K 



So at room temperature, we would expect both to rotate. 

It is interesting to note that if you put a bottle of hydrogen gas 
in liquid nitrogen (77 K), the molecules would hardly rotate at 
all! 
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Rotational heat capacity 



For most gases at room temperature, the energy levels are very 
closely spaced compared to k B T. The rotational partition 
function be computed using an integral: 

roc 0 

Z rot = / (2J+ l)exp[-h 2 J(J + l)/2Ik B T]dJ 
Jj=o 

The result is 

_ T 

z rot — q 

Note that this is a high temperature approximation, since we 
are assuming that k B T is much larger than the spacing between 
energy levels. 

The rotational energy at high temperature is then given by 

9 <9lnZ 

U = Nk B T 2 ^- = Nk B T 
and the rotational heat capacity by C ro t = Nk B . 
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From translation to rotation 



Imagine starting at a very low temperature for a diatomic gas, 
well below 0 rot . 

There would be little rotation, just translational motion. 

The heat capacity is the same as for monatomic gas: 

C = -Nk B = -nR. 

Increase the temperature gradually. When you come close to 
0 ro t, the some molecules start to rotate. 

This would take up additional energy. So the heat capacity 
increases. 

When the temperature gets much higher than 0 rot , the rotation 
would contribute an additional 

CVot = Nk B = R. 
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From rotation to vibration 



The spacing between energy levels for vibration tends to be 
much higher than those for rotation because of the strength of 
the bond. 



translation rotation vibration 
Following a similar reasoning to rotation, there would be a 
" characteristic temperature of vibration," # vib . This is also 
called the Einstein temperature. 

Below # vib , there would be little vibration. 

If the temperature is much higher than # vib , the vibration would 
contribute to the heat capacity. 
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Vibration 



The atoms in a diatomic molecule can vibrate because of the 
bod between them. We may imagine that it acts like a spring. 
The behaviour is approximately the same as the 1-D simple 
harmonic oscillator that we have seen previously. The energy 
levels are: 



E n = yn + - hv 

where v is the frequency. 

The same formula for the partition function, energy and heat 
capacity may be used. The heat capacity is given by 





where 6 — hu/k B . 
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Vibration 



The heat capacity is given by 

exp(fl/T) 
G-Nk B {-) (exp(0/T) _ 1)2 

where 0 = hvjk B . 



6 is defined to be # vib , the "characteristic temperature of 
vibration" : 

6> vjb = hu/k B . 

If the temperature is much higher than 0 vib , the vibration would 
contribute to the heat capacity an additional 

C - C vib = Nk B . 
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Temperature dependance 



Schematically, the heat capacity of a diatomic gas would vary 
with temperature in this way: 

C (R/mol) 



2.5 
1.5 

0 









/ 








J. 






1 


r 




Orot 


^vib 


T 



Experiments have shown that the heat capacities of real gases 
often follow this behaviour. 



There are some exceptions, like hydrogen 
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The heat capacity of hydrogen 



According to the theory, the heat capacity of hydrogen should 
follow curve A. 

2r t- -\ 
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Spin states of hydrogen 



The reason is because the protons in the nuclei of the molecule 
are fermions. 

The spin states of the protons affects the rotational quantum 
number J. This in turn affects the rotational energy. 

Protons are fermions. If the two protons are exchanged, the 
overall wavefunction must change sign. 

The overall wavefunction is a product of the rotation 
wavefunction and the spin wavefunction (as well as electronic, 
vibration, ...). 

When J is even, the rotation wavefunction remains the same 
when the protons are exchanged. So the spin wavefunction 
must change sign. This happens if the total spin is 0. 

When J is odd, the rotation wavefunction changes sign when 
the protons are exchanged. So the spin wavefunction must 
remain unchanged. This happens if the total spin is 1. 
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The equilibrium case 



The spin-0 molecule is called parahydrogen. 
The spin-1 molecule is called orthohydrogen. 

Assuming that the molecule can freely convert between the 
spin states, we can write down the partition function: 

Z rot = (2J+l)exp(-J(J+l)£ rot /^) 
J even 

+ 3 ^ (2J+l)exp(-J(J+l)0 rot /T) 
J odd 

where the factor of 3 comes from the degeneracy of the spin 
state 1. 

If we now compute the heat capacity ... 
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But they are not at equilibrium 



The reason for the discrepancy is that in hydrogen gas, the 
probability for a molecule to change from ortho-hydrogen to 
para-hydrogen (or vice-versa) in a collision with another 
molecule is extremely low. Molecules tend to stay in their spin 
states. 

If we take that into account, then we should treat the two 
components as different gases, each with their own rotational 
heat capacity 

If ortho-hydrogen has three-times greater concentration in a 
sample than para-hydrogen (because of the degeneracy of the 
spin-1 state), then the total heat capacity will be: 

_ 3 1 

Crot — -C'ortho + -Cpara 

This heat capacity gives the curve C in the previous figure... in 
agreement with experiment, at last! 
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What we have learnt so far. 



1. Density of states for particle in a box is 

2. The Maxwell-Boltzmann distribution breaks down when at 
low temperatures. 

3. "High temperature" means that k B T is much greater than 
the spacing between energy levels. 

4. The heat capacity for one mole of monatomic gas is 3R/2. 

5. For diatomic gas, rotation contributes another R to the heat 
capacity at high temperature. 



6. Vibration contributes a further R. 



7. The partition function formula for entropy is different for 
distinguishable and indistinguishable particles. 
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4.7 Exercises 
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Some exercises 



Exercise 1 



Show that the partition function of an ideal gas: 

roc AmrrA/ t 

Z= / ^^-(2ms) 1 / 2 exp(-s/^T)ds 
JO 

is given by: 

f27rmk B T\ 3 / 2 



Z = V 



[You are given that 



— oo n\(2a) Zn a 
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First, move the constant factors outside the integral: 

A-mnV 



Z = 



(2m) 1 / 2 J e 1/2 exQ(-e/k B T)de 



There is no power of 1/2 in the integration formula that is 
given. So we start with the substitution 



£ = X 2 



Then 

A-mnV 



~ (2m) 1 / 2 [°° xexp(-x 2 /k B T).2xdx 

JO 

(2m) 1 / 2 2 / x 2 exp(-x 2 /k B T)dx 
J 0 



h 3 

Next, set n to 1 in the given formula: 

•OO 2 _ a 2 x 2 27T 1 / 2 



— oo (2a)^a 

The lower of limit of this integral is — oo. To make this agree 
with the integral for the partition function Z ... 
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We change the lower limit to 0 and halve the right hand side, 
and the formula becomes: 

00 ^ cJjOC — 



'o (2a) 2 a 

Now comparing with the partition function again: 



Z = 



roc 

(2m) 1 / 2 2 / x 2 exp(-x 2 /k B T)dx 
JO 



the integrals become the same if we set 

2 1 
a = 



k B T 

Substituting this into the formula gives 

| o °° x 2 e-* 2 / k B T dx = ^(k B T) 3 / 2 
Substituting this into partition function intergral, 
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We get 



Simplifying, 



Z = -^(2m) 1 / 2 2^(^T) 3 / 2 



Z = v(«) 3/2 



we get the formula for the partition function 
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Exercise 2 



What is the root-mean-square speed of a helium atom at room 
temperature? 

[Given that the relative atomic mass of helium is 4. 
Atomic mass unit u is 1.6605 x 10 -27 kg. 
Boltzmann's constant k B is 1.3807 x 10~ 23 J K _1 . ] 



PHYS393/Hock 



84 



2009-2010 



First, we find the mean of v z . This is can be obtained by 



N Jo 

We have multiplied v 2 by the number particles at each energy 
interval de, and divided by the total number N. 

The energy of the atom in an ideal gas is just the kinetic 
energy, 



2 

If we rewrite the above integral for v 2 in terms of the energy, 
we get 




1 



2 



mv 




This integral is just the total energy U, i.e. 



2 — 



2 



V 



U. 



Nm 
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We have previously obtained the energy U from the partition 
function: 

U = -Nk B T. 
2 

So the mean square speed is 

-o 2 3 Ar7 _ 3k B T 

v 2 = Nk B T = — — 

Nm 2 m 

The root-mean-square speed is then 



v 2 = 



'3k B T 



m 
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We can now calculate the root-mean-square speed. The mass 
of the helium atom is 

Au = 4 x 1.6605 x 10~ 27 = 6.642 x 10~ 27 kg. 
Taking the room temperature as 298 K, the answer is 



3k B T 



m 



3 x 1.3807 x 10~ 23 x 298 
\ 6.642 x 10- 27 

1363 m/s 
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Exercise 3 



Given that the internuclear separation in the O2 molecule is 1.2 
A, calculate its characteristic temperature of rotation. 

[The relative atomic mass of the oxygen atom is 16. 
Atomic mass unit u is 1.6605 x 10 -27 kg. 
Planck's constant h is 1.0546 x 10 -34 J s. 
Boltzmann's constant k B is 1.3807 x 10~ 23 J K _1 . ] 
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The characteristic temperature of rotation is given by 
I is the moment of inertia, given by 

r would be the distance from the centre of the molecule to 
each atom, i.e. 

r = 1.2/2 = 0.6A, 

and the mass is 

m = 16u = 16 x 1.6605 x 10~ 27 = 2.657 x 10~ 26 kg. 
The O2 molecule has 2 atoms. So the moment of inertia is 

I = 2mr 2 = 2 x 2.657 x 10~ 26 x (0.6 x 10~ 10 ) 2 
= 1.913 x 10~ 46 kg m 2 . 
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The characteristic temperature of rotation is given by 

h 2 



#rot — 



Ik B 

(1.0546 x 10~ 34 ) 2 

1.913 x IO- 46 x 1.3807 x 10~ 23 
4.21K. 
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Exercise 4 



The Einstein temperature (the characteristic temperature of 
vibration, # vib ) of O2 is 2200 K. Calculate the percentage 
contribution of the vibrations to its heat capacity at room 
temperature. 
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The energy levels of the 1-D simple harmonic oscillator are 
given by 

E n = (n+l/2)hv. 

where v is the frequency of vibration. The Einstein 
temperature is defined by 

#vib = hv/k B , 

As room temperature is much lower than this, we may expect 
most of the molecules would at at the lowest level E 0 . Perhaps 
a very small number would be excited to the next level E\. 

We can use this idea to estimate the vibrational heat capacity 
at room temperature. 
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So we would be interested into only 2 levels, E 0 and E 1 . The 
spacing in between is hv. 

The probability of a particle at one level is proportional to the 
Boltzmann factor. Since most of the particles are at Eq, the 
probability of finding a particle there is close to 1. So the 
probability of finding a particle at E\ would be given by the 
Boltzmann factor exp(-hv/k B T). 

This means that the number of particles in: 

Ei is ni — N exp(—hv/k B T), and 
E 0 \s n 0 = N - Nexp(-hv/k B T). 
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Because exp(—hv/k B T) would be very small, the total is still 
close to N. The total energy is then 

U = n 0 E 0 + n 1 E 1 

1 3 

= -hv{N - Nexp(-hv/k B T)) + -hvN exp(-hv/k B T) 

1 

= -hvN{l + 2exp{-hv/k B T)) 
The heat capacity is then 

c = dU 



dT 



= hvN(hv/k B T z )exp(-hv/k B T)) 
= Nk B (hu/k B T) 2 exp(-hv/k B T)) 
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In terms of the Einstein temperature, 

#vib = hu/k B , 

the heat capacity is 

C = Nk B (0 ylb /T) 2 exp(-0 ylb /T)) 

This would be valid at temperatures much lower than # vib , such 
as room temperature. At 298 K, the vibration heat capacity is 
then given by 

C vib = Wc B (2200/298) 2 exp(-2200/298) = 0.0339Nk B 

At room temperature, the heat capacity mainly consists of 
translational and rotational contributions, which gives 

C = ^Nk B + Nk B = ^Nk B 

Therefore the percentage contribution of the vibrations to the 
heat capacity at room temperature is 

0.0339^- = 13.6 
2 
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